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Abstract. Cross sections for the rotational (de)excitation of CO by ground state para- and ortho-H 2 are obtained using quantum 
scattering calculations for collision energies between 1 and 520 cirT 1 . A new CO-H2 potential energy surface is employed and 
its quality is assessed by comparison with explicitly correlated CCSD(T)-R12 calculations. Rate constants for rotational levels 
of CO up to 5 and temperatures in the range 5-70 K are deduced. Fhe new potential is found to have a strong influence on the 
resonance structure of the cross sectio ns at very low collision energies. As a result, the present rates at 10 K differ by up to 50% 
with those obtained bv lFlowedl200ll) on a previous, less accurate, potential energy surface. 
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1. Introduction 

Since its discovery i n interstellar space more than 30 years ago 
ftVilson et al] fl970). carbon monoxide (CO) has been exten- 
sively observed both in our own and in other galaxies. CO 
is indeed the second most abundant molecule in the Universe 
after H2 and, thanks to its dipole moment, it is the principal 
molecule used to map the molecular gas in galactic and extra- 
galactic sources. CO is also currently the only mol ecule to have 
been observed at redshifts as high as z — 6.42 JWalter et alJ 
I2003I) . Collisional excitation of rotational levels of CO occurs 
in a great variety of physical conditions and emission from lev- 
els with rotational quantum numbers J up to 39 have been id en- 
tified in circumstellar environments ( Cernic haro et all 1 996 ) . 

Hydrogen molecules are generally the most abundant col- 
liding partners in the interstellar medium, although collisions 
with H, He and electrons can also play important roles, for ex- 
ample in diffuse interstellar clouds. Rate coefficients for colli- 
sions with H2 based on the potential energy surface (PES) of 
Jankowski & Sz alewicz (199 8, hereaf ter JS98) ha ve been cal- 
culated recently dMengeletalJ200lHFlowerl200ll) . These rates 
were found to be in reasonable agreement (typically within a 
factor of 2 at 10 K) with those comp uted on older and less ac- 
curate PES (e.g. lSchinke et alJ dl985t) '). H owever, in con trast to 
previous works. iMmgeTetaTI l l200ll) and lFloweil J200ll) found 
the inelastic rates with AJ = 1 larger than those with AJ — 2 
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for para-H2. This result reflects the crucial importance of the 
PES, whose inaccuracies are one of the main sources of er- 
ror in collision al rate ca lculations, especi ally at low tempera- 
tures (see, e.g., Balakrishnan et al. (2002)). The JS98 surface 
has been checked ag a inst ex periment by a couple of authors. 
iGottfried & McB and (120011) . in particular, computed second 
virial coefficients below 300 K and showed that the attractive 
well of the PES was slightly too deep. A modification which 
consists of multiplying the PES by a factor of 0.93 was there- 
fore suggested b y these authors. This modification was applied 
bv lMengel et~al1J200lb and it was indeed found to produce bet- 
ter agreement with measurements of pressure broadening cross 
sections. On the other hand, the original PES was found to give 
a very good agreement with experimental stat e-to-state rota- 
tiona l cross sections measured near 1000 crrr 1 llAntonova et al.l 
2000). At these energies, however, onl y the short-ran ge repul- 
sive region of the PES is probed. Finally, Fl ower! d200ll) showed 
that the scaling of 0.93 has only minor influence on the rota- 
tional rates and he employed the original PES. 

Very recently, a new P ES has been computed by 
Llankowski & Szalewicz! d2005l) (hereafter JS05). In contrast to 
the JS98 surface, which was calculated within the symmetry- 
adapted perturbation theory (SAPT), the JS05 surface has been 
obtained using the supermolecular approximation based on the 
CCSD(T) method using correlation-consistent basis sets com- 
plemented by bonding functions and basis set extrapolation. It 
turned that the major inaccuracy of JS98 was due to the use of 
a correction for the effects of the third and higher orders ex- 
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tracted in an approximate way from supermolecular SCF cal- 
culations. It recently became clear that this correction should 
not be used for nonpolar or nearly nonpolar monomers. If this 
term is removed from JS98, the agreement of the re sulting 
surface with JS05 is excellent. Jank owski & Szalewiczl J2005) 
have employed their new PES to calculate rovibrational energy 
levels of the para-H2-CO complex and they obtained a very 
good agreement (better than 0.1 cirT 1 comp ared to 1 cm" 1 fo r 
the JS98 surface) with the measurements by M cKellaJ (H"998). 
The computed values of the second virial coefficients also agree 
ve ry well with the experimental on es. Full details can be found 
in I.Tankowski & Szalewiczl ll2005h . In the present paper, we 
present new rate constants for the rotational excitation of CO 
by H2 based on the JS05 surface. Our results are restricted to 
low temperatures (T<70 K) because the new PES calculations 
are found to affect the rotational cross sections only at low col- 
lision energy (E rr ,n < 60 cm" 1 ). As shown below, differences 
with the results of Flower (2001) range from a few percent to 
50 %. Such effects, if not dramatic, are significant in view of 
the special importance of the CO-H2 system in astrophysics. 
Furthermore, the present results are of particular interest for 
comparisons with future measurements planned in Rennes by 
Sims and collaborators (Sim sl20o'5l) . 



2. Method 

2.1. Potential Energy Surface 

Following the conventions of JS98, the CO-H2 "rigid-body" 
PES is described by the following four Jacobi coordinates: 
three angles 0\, 82 and <p and the distance R denoting the inter- 
molecular centre of mass separation. The angles 9\ and 62 are 
the tilt angles of H2 and CO with respect to the intermolecu- 
lar axis, and (p is the dihedral angle. In the linear geometry, the 
oxygen atom points toward H2 when #2=0. The JS05 surface 
has been obtained by averaging over the intramolecular ground 
state vibration of H2, while CO was treated as rigid at its vibra- 
tional ground-state averaged geometry. The JS05 surface has 
a global depth of -93.049 cm" 1 for the linear geometry with 
the C atom pointing toward the H2 molecule (0\ = 0, 02=180°, 
= 0) and for a centre of mass separation R m \ n =7.92 bohrs. 
The (unsealed) JS98 and JS05 surfaces are presented in Fig. 1 
and compared to independent CCSD(T)-R12 calculations us- 
ing ground-state averaged geometries for CO and H2. 

The latter calculations can be considered as high accuracy 
state-of-the-art ab initio calculations. Indeed the CCSD(T)- 
R12 approach 1 has been shown to produce intermolecular en- 
ergies with n ear spectroscopic accuracy without recourse to 
extrapolation dHalkier et al.lll999t iMiiller & Kutzelniggfe oOO ; 
iFaure et alJ2005l) . In order to save some computing time and to 
improve the numerical stability of the R12 operator, the core or- 
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1 The R12 coupled cluster theory with singles, doubles and pertur- 
bative triples properly describes the electron-electron cusp by includ- 
ing explicitly the int er-electronic coordin ates into the wav efunction 
iNoaa & Kutzelniae 1 1994), see also lNoea & ValirorJ EooH) for a re- 
view.) 
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Fig. 1. H2-CO interaction energy as a function of R for three 
coplanar configurations: (6> 1 ,6» 2 )= (90°, 90°), (0°, 0°) and (0°, 
180°). The dashed and solid lines denote the JS98 and JS05 
calculations, respectively. The filled circles correspond to in- 
dependent CCSD(T)-R12 calculations (see text). 

bitals were frozen 2 in the CCSD(T)-R12 calculations. We em- 
ployed R12-suited basis sets including up to g funct ions for C 
and O , and d functions for H (as developed in lNoga & ValirorJ 
d2003l) : iKedzuch et al. (2005)). The R12 results are assumed 
to approach the basis set limit within 1-2 cm -1 . In particular 
the inclusion of a large set of bonding functions (corresponding 
to the R12-suited H set) affects the interaction by as 0.5 cm" 1 
only at 7 bohrs separation, while the counterpoise corrections 
remain very small, generally below 1 cm" 1 . 

Here, we can observe that the JS05 surface agree very 
well with the CCSD(T)-R12 results while the JS98 surface 
overestimates the depth of the CO-H? potential well by 
about 10-20 %, as anticipated bv lGottfried & McBanelll200ll) . 
However, the JS05 surface lies slightly below the CCSD(T)- 
R12 one. This agrees with the observation in JS05, based on 
extrapolations involving bases up to quintuple zeta quality, that 
the potential is about 1-2 cm" 1 too deep in the region of the 
global minimum. On the other hand, the comparison of the 
virial coefficients from JS05 with experiment suggests that the 
exact potential is still deeper, probably due to the effects re- 
sulting from theory levels beyond CCSD(T). We also checked 
that the R12 values are recovered within 1 cm" 1 by a basis set 
extrapolation involving larger (triple and quadr uple zeta) sets 
and th e popular bonding functions proposed by Will iams et alJ 
dl995l) . These comparisons provides a further check of the 
quality of the new CO-H2 PES and also confirm the suspected 
inaccuracies of the JS98 surface. 

In order to interface the JS98 and JS05 surfaces 
with the scattering co de, t he Fortran routines of 
Llankowski & Szalewiczl Jl998h : I.Tankowski & Szalewiczl 
ll2005h were employed to generate a grid of 57,000 geometries. 
These grid points were distributed in 19 fixed intermolecular 
distances R (in the range 4-15 bohrs) via random sampling 



- This approximation is of little consequence as core correlation ef- 
fects are minor for the CO-H 2 interaction (about 0.3 cm" 1 at 7 bohrs 
separation). 
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for the angular coordinates. At each intermolecular distance, 
the interaction potential V was least-squares fitted over the 
following angular expansion: 



100 



V(R,e u 02,4>) = YjV lM (R)s llhl (eue 2 , 
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where the basis functions s^y are produ cts of spher ical har- 
monics and are explicited in Eq. (A9) of Green ( 1975). These 
functions form an orthonormal basis set, and the expansion co- 
efficients Vi.ui can thus be written as 3 



CD 
CO 
CO 
CO 

o 

o 



v hh i(R) 



-iU: 

Je, Jfli Jd> 



V(R,6 1 ,e 2 ,(p)sf l ,(01,02,0) sinM#i sint9 2 <ie 




100 



The latter expression remains valid for any truncated expan- 
sion. 

We selected an initial angular expansion including all 
anisotropies up to Zi = 10 for CO and ^2=6 for H2, resulting in 
142 Si t i 2 i functions. The accuracy of the angular expansion was 
monitored using a self-consistent Monte Carlo error estimator 
jRist et alJ2 005). which also permitted to select the most perti- 
nent angular terms. The resulting set was composed of 60 basis 
functions involving anisotropies up to l\ -1 and l 2 — 4. A cu- 
bic spline interpolation was finally employed over the whole in- 
termolecular distance range and was smoothly connected with 
standard extrapolations to provide continuous expansion coef- 
ficients suitable for scattering calculations. The final accuracy 
of this 60 term expansion was found to be better than 1 cm 
in the global minimum region of the PES while the individ- 
ual expansion coefficients were converged within 0.01 cm -1 . It 
should be noted that the expansion does not reproduce exactly 
the values of the potential energ y at the grid p oints 4 , in con- 
trast to the approach followed by Flower (2001) who used the 
same number (25) of angular geometries and basis functions. 
The main advantage of our approach is to guarantee the con- 
vergence of the individual expansion coefficients, in particular 
in the highly anisotropic short-range part (R < 6 bohrs) of the 
interaction. 



2.2. Scattering calculations 

The molscat code JHutson & Greenl ll994t) was employed to 
perform the rigid-rotor scattering calculations reported below. 
All calculations were made using the rigid rotor approximation, 
with rotational constants Bco = 1.9225 cm -1 andfin, =60.853 
cnT 1 . For para-H2 (j 2 = 0), calculations were carried out at 
total energies (collision plus CO rotation) ranging from 4 to 
520 cm -1 . Full close-coupling calculations were performed be- 
tween 4 and 160 cirT 1 , with steps of 0.2 cirT 1 in the resonance 
range 4-120 cirT 1 , and a much coarser step of 5 cirT 1 between 
120 and 160 cm . Between 160 and 520 cm , the coupled- 
state approximation ( McGuir e & Kouril 19741) was shown to re- 
produce full close-coupling calculations within a few percent 

3 ★ denot es the complex conjugate. Note that the angular con- 
ventions of iGree nl i 1975§> are slightly different from those of 



ventions of Green 1 1975) ar 
Ijankowski & SzalewiczH i998). 
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4 The residual differences indicate the cumulative effects of the 
higher expansion terms. 



Fig. 2. Inelastic deexcitation cross sections for CO in collisions 
with para-H2 as a function of total energy. Both curves are the 
results of our calculations with the unsealed JS98 PES (dashed 
curves), and the JS05 PES (solid curves). 

and was thus employed. The energy step was 20 crrr 1 between 
160 and 340 crrr 1 , and 30 cirr 1 at higher energies. For ortho- 
H2, the same collision energy grid was used, resulting in to- 
tal energies in the range 126-64 2 cm' 1 . The propagator used 
was the log-derivative propagator ( Alexander & Manolopoulos 
1987). Parameters of the integrator were tested and adjusted to 
ensure a typical precision of 1-2 %. All the calculations with 
para-H2 used a j 2 = 0, 2 basis, while the ortho-H2 calculations 
used only a j 2 = 1 basis. Tests of this basis were made at the 
highest energy, where inclusion of the j 2 = 3 (closed) state was 
found to affect cross sections by less than one percent. At least 
three closed channels were included at each energy for the CO 
rotational basis. At the highest energy, we had thus a j\ < 18 
basis. Transitions among all levels up to j\ - 5 were computed. 
Finally, rate constants were calculated for temperatures ranging 
from 5 to 70 K by integrating the cross sections over essentially 
the whole collision energy range spanned by the corresponding 
Maxwell-Boltzmann distributions. The energy grid was cho- 
sen so that the highest collision for all transitions be about ten 
times larger than the highest temperature. The accuracy of the 
present rates is expected to lie between 5 and 10 % at 70K, and 
probably better at lower temperatures. 

3. Collisional cross sections and rates 

Figure |3 clearly illustrates the fact that the resonance pattern, 
which is here of 'shape' rather than 'Feshbach' type as dis- 
cussed by iFlowerl d200ll) . strongly depends on the interaction 
potential. We also note that the very fine grid used at the lowest 
energies was necessary for a proper description of the reso- 
nances. It can be observed that the cross sections converge at 
higher energies (above the resonance regime), as the collision 
energy becomes large compared to the potential energy well 
depth. As a result, ra te constants are expected to differ signif- 
icantly from those of iFlowerl d200lh at very low temperatures 
only. 

We present in Tables [T] and |3] our rate constants at 10, 30 
and 70 K, and those of Flower (2001) for comparison. We note 
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that especially at 10 K, rates can vary up to 50%. These differ- 
ences reflect the influence of three independent parameters: (i) 
the high accuracy of the new ab initio potential energy calcula- 
tions, (ii) the use of a well converged expansion of the PES over 
the angular basis functions and (iii) the careful description of 
resonances. It is however not clear how to discriminate among 
the three effec ts. At 70K (and beyond), our rates and those of 
iFlowerl J200lh agree within typically 10-20 %, as anticipated 
from the convergence of the cross sections at high energy. We 
therefor e conclude tha t the differences between our results and 
those of Flowei] J200lh are astrophysically relevant for temper- 
atures below 70 K. Above this value, the inaccuracies of the 
JS98 PES have on ly a minor infl uen ce on the rotational rates 
and the results of lFloweJ fcOOll) and lMeneel et alJ J200ll) are 
reliable. 

The para-H2 rates presented in the 30 K panel in Table \l\ 
illustrate the respective contribution of the surface and of the 
subsequent angular expansion and collisional treatment. Data 
in columns 1 and 2 are based upon the same JS98 surface 
and the differences are ascribable to the better convergence 
of the angular expansion and to the finer description of the 
resonances. Conversely, data in columns 2 and 3 are based 
upon the same angular expansion and collisional treatment, 
and the differences only reflect the improvement in the surface. 
Depending on initial and final rotation states, various situations 
may occur. In some cases there is a compensation of errors for 
the surface (as for 2-1, 3-0, 3-2), or for the expansion and col- 
lisional treatment (as for 5-0). Collisional and surface errors 
may also partially cancel out, as for the 1-0 rate, or add up as 
for the 2-0 rate. While the detailed behaviour seems complex 
and irregular, in average the errors introduced by the older JS98 
surface are of the same order of magnitude as the errors intro- 
duced by the limited a ngular expansion and energy mesh in the 
older work by Flower ( 2001 ). 

It is interesting to note that for collisions with para-H2, in- 
elastic rates with Aj\ = 1 are larger than those with Aj\ -2 
while the reverse applies for ortho-H^. The propensity to fa- 
vo r even A n over odd A/i h as been explained semiclassically 
bv lMcCur dv & Miller! i 1977b in terms of an interference effect 
related to the even anisotropy of the PES. The reverse propen- 
sity can also occur if the odd anisotropy is suffici ently large, 
as ob served experimentally in the case of CO-He JCartv et alJ 
120041) . Here, the propensity depends on the quantum state of 
the projectile because in the para case (y'2 = 0), some terms 
of the PES expansion (those associated with the quadrupole of 
H2) vanish identically. Therefore, the even anisotropy of the 
PES is larger with ortho-H2 than with para-H2 and Aji -2 are 
favored. Note, however, that these propensity rules depend on 
j\, Aj\ and the temperature. An other inte resting point is that 
the differences with the results of lFloweJ fooOli are generally 
smaller with ortho-H2 than with para-H2. This again reflects 
the different anisotropies of the PES for para and ortho-H2: in 
the ortho case, cross sections are larger and scattering calcu- 
lations are less sensitive to small changes in the PES. Similar 
effe cts have also been ob served recently in the case of H2O-H2 
jDubernet&etalJ2004. 

It is also of some interest to co mpare the present re- 
sults with the CO-He rates obtained bv lCecchi-Pestellini et alJ 



(2002|). It is generally assumed that rates with para-H2 (y'2 = 0) 
should be about 50% larger than He rates owing to the smaller 
collisional reduced mass and the larger size of H2. We have 
checked that the present H2 (72 = 0) rates are in fact within 
a factor of 1—3 of the He rates. Thus, as already shown by 
Philli£S£LalJ|l996) in the case of H2O, rates for excitation by 
H2 are not adequately represented by scaled rates for excitation 
by He. 

Finally, for use in astrophysical modelling, the temperature 
dependence of the above tra nsition rates k(T) have bee n fitted 
by the analytic form used by Balakrish nan et al] |l999): 

4 

logmen = 2] ( 3 ) 

«=o 

where the x factor was changed from r~'/ 3 to r~'/ 6 to achieve 
a better fitting precision. The quadratic error on this fit is lower 
than 0.1% for all transitions. The coefficients of this fit are pro- 
vided as online material. For those who plan to use them, we 
emphasize that these fits are only valid in the temperature range 
5 < T < 70 K. Finally, excitation rates are not given in this pa- 
per, but can be obtained by application of the detailed balance. 

4. Conclusion 

Cross sections for the rotational (de)excitation of CO by ground 
state para- and ortho-H2 have been computed using quan- 
tum scattering calculations for collision energies in the range 
1-520 cm . A new, highly accurate, CO-H2 potential en- 
ergy surface has been employed and it has been shown to 
strongly influence the resonance structure of the cross sections 
in the very low collision energy regime (E co n < 60 cm -1 ). 
Conversely, at higher energies, the effect of the new poten- 
tial was found to be only minor. As a result, the pres ent rate 
consta nts are found to differ significantly from those of Fl owed 
J200ll) . obtained on a previous and less accurate potential, only 
for temperatures lower than 70 K. Transitions among all levels 
up to ji - 5 only were computed as higher states are generally 
not populated at the low temperatures investigated here. 

This work illustrates the relationship between the accuracy 
of the potential energy surface and the accuracy of the corre- 
sponding inelastic cross sections for a simple system of astro- 
physical relevance. Even at low temperatures, moderate inac- 
curacies in the surface (in the 10-20% range) result in semi- 
quantitative inelastic rates with typical errors in the 20-50% 
range with no dramatic amplification of errors. However, this 
optimistic conclusion presents also a harder counterpart, as we 
show that in order to obtain accurate inelastic rates one has to 
satisfy all three conditions, i.e., use an accurate surface, a prop- 
erly converged angular expansion, and a detailed description of 
the cross-section resonances. Consequently, the improvement 
of the accuracy of inelastic rates, especially at higher tempera- 
tures when the number of channels is large, represents a consid- 
erable computational effort. In this respect, the work by Fl oweJ 
(2001) represented an excellent compromise to achieve a typi- 
cal 20% accuracy below 100 K. 

This work extends the initial objectives of the "Molecular 
Universe" European FP6 network (2005-2008) and participates 
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to the current international efforts to improve the description of 
the underlying microscopic processes in preparation of the next 
generation of molecular observatories. 

Acknowledgements. CCSD(T)-R12 calculations were performed on 
the IDRIS and CINES French national computing centers (projects 
no. 051141 and x2005 04 20820). MOLSCAT calculations were per- 
formed on an experimental cluster using a subset of the computer 
grid tools (under project Cigri of the "Action Incitative GRID") with 
the valuable help from F. Roch and N. Capit. This research was 
supported by the CNRS national program "Physique et Chimie du 
Milieu Interstellaire" and the "Centre National d'Etudes Spatiales". 
MW was supported by the Ministere de l'Enseignement Superieur et 
de la Recherche. KS acknowledges the NSF CHE-0239611 grant. 



Phillips, T. R., Maluendes, S., & Green, S. 1996, ApJS, 107, 
467 

Rist, C, Faure, A., & Valiron, R 2005, in preparation 
Schinke, R., Engel, V., Buck, U., Meyer, H., & Diercksen, 

G. H. F. 1985, ApJ, 299, 939 
Sims, I. R. 2005, private comunication 

Walter, F, Bertoldi, F, Carilli, C., et al. 2003, Nature, 424, 406 
Williams, H. L., Mas, E. M., Szalewicz, K., & Jeziorski, B. 

1995, J. Chem. Phys., 103, 7374 
Wilson, R. W., Jefferts, K. B., & Penzias, A. A. 1970, ApJ, 161, 

L43+ 



References 

Alexander, M. H. & Manolopoulos, D. E. 1987, J. Chem. Phys., 
86, 2044 

Antonova, S., Tsakotellis, A. P., Lin, A., & McBane, G. C. 

2000, J. Chem. Phys., 1 12, 554 
Balakrishnan, N., Forrey, R. C, & Dalgarno, A. 1999, ApJ, 

514,520 

Balakrishnan, N., Yan, M., & Dalgarno, A. 2002, ApJ, 568, 443 
Carty, D., Goddard, A., Sims, I. R., & Smith, I. W. M. 2004, 

J. Chem. Phys., 121,4671 
Cecchi-Pestellini, C, Bodo, E., Balakrishnan, N., & Dalgarno, 

A. 2002, ApJ, 571, 1015 
Cernicharo, J., Barlow, M. J., Gonzalez-Alfonso, E., et al. 

1996, A&A, 315, L201 
Dubernet, M. L. & et al. 2005, in preparation 
Faure, A., Valiron, P., Wernli, M., et al. 2005, J. Chem. Phys., 

122, 221102 

Flower, D. R. 2001, J. Phys. B: At. Opt. Mol. Phys., 34, 2731 
Gottfried, J. & McBane, J. C. 2001, J. Chem. Phys., 1 12, 4417 
Green, S. 1975, J. Chem. Phys., 62, 2271 
Halkier, A., Klopper, W., Helgaker, T., J0rgensen, P., & Taylor, 

P. R. 1999, J. Chem. Phys., 1 1 1, 9157 
Hutson, J. M. & Green, S. 1994, MOLSCAT computer 

code, version 14 (1994) (distributed by Collaborative 

Computational Project No. 6 of the Engineering and 

Physical Sciences Re search Council UK) 
Jankowski, P. & Szalewicz, K. 1998, J. Chem. Phys., 108, 3554 
Jankowski, P. & Szalewicz, K. 2005, The Journal of Chemical 

Physics, 123, 104301 
Kedzuch, S., Noga, J., & Valiron, P. 2005, Mol. Phys., 103, 999 
McCurdy, C. W. & Miller, W. H. 1977, J. Chem. Phys., 67, 463 
McGuire, P. & Kouri, D. J. 1974, J. Chem. Phys., 60, 2488 
McKellar, A. R. W. 1998, J. Chem. Phys., 108, 1811 
Mengel, M., de Lucia, F. C, & Herbst, E. 2001, Canadian 

Journal of Physics, 79, 589 
Miiller, H. & Kutzelnigg, W 2000, Phys. Chem. Chem. Phys, 

2, 2061 

Noga, J. & Kutzelnigg, W. 1994, J. Chem. Phys., 101, 7738 
Noga, J. & Valiron, P. 2002, in Computational Chemistry: 

Reviews of Current Trends, Vol. 7 (J.Leszczynski ed., World 

Scientific, Singapore), 131 
Noga, J. & Valiron, P. 2003, Collect. Czech. Chem. Commun, 

68, 340 



6 



M. Wernli et al.: improved low-temperature CO-H 2 rotational rate constants 



Table 1. Rotational deexcitation rates of CO in collision with 
para-H2 (72 = 0) at three temparatures. I and F are the CO ini- 
tial and final rotation states, respectively. All rates are given in 
units of 10~ 10 cm 3 /s. Third column recalls the rates of Flower 
obtained with the unsealed JS98 surface. Fourth col- 
umn presents our rates obtained with the JS05 surface. In the 
30 K panel, the rates obtained with our expansion of the JS98 
surface are also presented. 
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2 


0.153 


0.125 








4 


3 


0.888 


0.822 








5 


3 


0.656 


0.759 








5 


4 


0.925 


0.662 








I 


F 


Flower 


This work 












Temp.: 70 K 










1 





0.336 


0.344 








2 





0.300 


0.320 








3 





0.098 


0.083 








4 





0.035 


0.044 








5 





0.020 


0.024 








2 


1 


0.636 


0.602 








3 


1 


0.486 


0.509 








4 


1 


0.177 


0.161 








5 


1 


0.061 


0.080 








3 


2 


0.681 


0.659 








4 


2 


0.567 


0.591 








5 


2 


0.221 


0.205 








4 


3 


0.678 


0.689 








5 


3 


0.631 


0.654 








5 


4 


0.748 


0.705 




I 


F 


Flower 


This work JS98 This work JS05 






Temp.: 30 K 






1 








.300 


0.288 




0.327 


2 








.269 


0.297 




0.315 


2 


1 





661 


0.613 




0.612 


3 








.072 


0.059 




0.060 


3 


1 





.456 


0.492 




0.509 


3 


2 





.695 


0.677 




0.670 


4 








.026 


0.029 




0.037 


4 


1 





.130 


0.110 




0.118 


4 


2 





.555 


0.604 




0.612 


4 


3 





.747 


0.716 




0.712 


5 








.010 


0.010 




0.014 


5 


1 





.047 


0.054 




0.068 


5 


2 





.175 


0.149 




0.160 


5 


3 





.648 


0.712 




0.705 


5 


4 





.816 


0.740 




0.713 



Table 2. Rotational deexcitation rates of CO in collision with 
ortho-H2 (j'2 = 1)- The conventions are the same as for the para 
case. 



I 


F 


Flower 


This work 






Temp.: 10 K 




1 





0.388 


0.379 


2 





0.596 


0.559 


3 





0.105 


0.071 


4 





0.083 


0.061 
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0.015 
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z 


1 

1 


U. / JVJ 


719 


1. 

j 


1 

1 


f) 84? 


O 831 


A 
L + 


1 

1 


f) 1 SO 


O 1 94 


J 


1 

1 


O 1 94 


n 1 


"2 
J 


7 

J. 


n soi 


n 74.8 


A 
~t 


9 

Z 


n 86? 

u.ouz 


f) 804. 

U.07t 


D 


9 
Z 


U.loo 


n 1 


A 
t 


J 


U.oZ / 


n 786 

u. / OU 


J 


J 


O 8Q6 


fl Q79 


< 
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871 
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I 


F 


Flower 


This work 






Temp.: 70 K 




1 





0.402 


0.350 


2 





0.629 


0.468 


3 





0.114 


0.090 


4 





0.079 


0.063 


5 





0.027 


0.028 


2 


1 


0.684 


0.620 


3 


1 


0.858 


0.728 


4 


1 


0.199 


0.174 


5 


1 


0.124 


0.111 


3 


2 


0.732 


0.686 


4 


2 


0.871 


0.808 


5 


2 


0.236 


0.216 


4 


3 


0.757 


0.711 


5 


3 


0.880 


0.844 


5 


4 


0.779 


0.719 
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I 


F 


a 


«i 


«2 


a 3 


a 4 


1 





-10.413 


4.162 


-20.599 


32.563 


-16.889 


2 





-9.202 


-7.997 


18.452 


-18.717 


6.910 


2 


1 


-0.379 


-59.429 


131.662 


-127.161 


45.605 


3 





-3.319 


-35.329 


53.745 


-29.641 


2.471 


3 


1 


-7.042 


-20.755 


48.788 


-49.875 


18.607 


3 


2 


-1.976 


-48.986 


107.139 


-102.134 


36.288 


4 





-8.004 


-16.164 


28.359 


-21.986 


6.026 


4 


1 


1.979 


-69.801 


140.821 


-125.047 


40.556 


4 


2 


-5.102 


-34.014 


82.717 


-87.211 


33.665 


4 


3 


-1.891 


-49.861 


109.419 


-103.434 


35.726 


5 





7.564 


-106.151 


220.187 


-205.770 


72.462 


5 


1 


-5.209 


-32.541 


67.718 


-63.996 


23.105 


5 


2 


4.736 


-91.653 


203.779 


-201.295 


73.681 


5 


3 


-5.607 


-31.380 


78.077 


-83.702 


32.930 


5 


4 


-0.510 


-60.401 


137.793 


-134.344 


46.449 



Table .1. Fitting coefficients of CO-para-H2 rates, following 
formula |3 / and F are the initial and final rotation states, re- 
spectively. The rates thus obtained are in cm 3 /s. 



/ 


F 


flo 


«i 


ai 


«3 


«4 


1 





-7.351 


-11.062 


5.033 


16.394 


-14.176 


2 





0.331 


-68.556 


162.123 


-167.243 


63.863 


2 


1 


4.370 


-90.464 


207.213 


-207.995 


77.557 


3 





4.128 


-88.472 


191.123 


-181.945 


64.170 


3 


1 


-1.466 


-57.252 


138.316 


-144.752 


55.539 


3 


2 


-0.194 


-61.716 


141.443 


-142.699 


53.788 


4 





-6.118 


-31.232 


70.967 


-70.685 


25.956 


4 


1 


5.330 


-96.359 


215.047 


-212.299 


77.662 


4 


2 


-2.940 


-49.481 


124.615 


-135.311 


53.609 


4 


3 


-1.454 


-53.713 


122.102 


-121.151 


44.472 


5 





5.241 


-96.162 


206.655 


-199.255 


72.234 


5 


1 


-6.365 


-27.561 


61.697 


-61.299 


22.802 


5 


2 


6.689 


-107.753 


249.875 


-256.409 


97.478 


5 


3 


-3.812 


-45.669 


119.845 


-134.502 


54.860 


5 


4 


1.082 


-71.954 


169.486 


-173.174 


64.402 



Table .2. Fitting coefficients of CO-ortho-H2 rates, following 
formula / and F are the initial and final rotation states, re- 
spectively. The rates thus obtained are in cm 3 /s. 



